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Abstract 



%4 ' We compute analytically the effects of energy conservation on the self-similar structure 

of parton correlations in QCD jets. The calculations are performed both in the constant and 
running coupling cases. It is shown that the corrections are phenomenologically sizeable. 
On a theoretical ground, energy conservation constraints preserve the scaling properties of 
correlations in QCD jets beyond the leading log approximation. 
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0\ ■ 1 Introduction 

m . 

The structure of parton multiplicity correlations within QCD jets and in particular their self- 
similar properties have been recently studied 0, g f§. The main result is twofold : 

& ' • The structure of a QCD jet is multifractal, in the sense that the multiplicity density of 

gluons does not occupy uniformally the available phase-space when it is computed in the 



Qh. Double-Leading-Log approximation (DLA) of the perturbative QCD expansion. More 

precisely, if one measures the inhomogeneities of the multiplicity distribution using the 
scaled factorial moments of order q, one writes : 

>< [ F q (A) ee (^-1) (n - q + l)) A ^ A ^ m ^ Vq/d) (1) 

bj , \ n /A 

where n is the multiplicity of partons registered in a small phase-space interval A. T> q 
is called the fractal dimension of rank q of the density of gluons while d is the overall 
dimension of the phase space under consideration (in practice, d = 2 for the solid angle, 
and d = 1 if one has integrated over, say the azimutal angle, by keeping fixed the opening 
angle with respect to the jet axis). The multiplicity distribution is uniform if T> q = d, 
while it is fractal if it is smaller. Then multifractality is for a q— dependent dimension. 

In the DLA, assuming a fixed QCD coupling constant a s , the result is the following: 

V q = 70 — (2) 
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where 70, considered to be small, is given by : 

7o = 4iV c gj. ; N c = 3 (3) 

In the running coupling constant scheme, the jet structure is modified by scaling violation 
effects. It leads to a multifractal dimension slowly varying with the angular variables 
of observation. Furthermore, one observes^, at small angles V g — ► d, i.e. a dynamical 
saturation of multifractality. The calculation gives : 

V q {z) ~ 7 o (a,(£7flo)) ^ l - ^ _ = (4) 



where the scaling variable z is defined as: 



log^o/A 



(5) 



Note that the fractal dimension is no more constant and depends on both the angular 
direction(#o) and angular aperture {&) of the observation window. In practice, the fractal 
dimensions have a tendency to increase with z and 70 and to reach the saturation point 
where they become equal to the full dimension d. This saturation effect comes from the 
increase of the coupling constant at larger distances and signals the onset of the non- 
perturbative regime of hadronization. 

We want to reconsider these results by going beyond the approximations made in Refs.|l], |2|, 4|. 
Some results exist which include Next-Leading-Log corrections to the anomalous dimensions[l|, 
but they do not include the energy-conservation constraints. Indeed, it is well known that per- 
tubative QCD resummation in the leading log approximation (DLA) predict too strong global 
multiplicity moments, at least when compared to the hadronic final state observed in experimen- 
tal data on jets||. In fact, as pointed out in Ref.|6|, the energy conservation (EC) constraint 
at the triple parton vertices may explain the damping of multiparticle moments observed in the 
data. While the EC correction is perturbatively of higher order {a s as compared to and 
gives a rather moderate effect on the mean multiplicity, the calculation of multiplicity moments 
of higher rank q gives a contribution of order q 2 o?J 2 which become important already for q = 2. 
As we shall demonstrate, this effect is even more important for local correlations in phase-space. 

The main goal of our paper is the analytical calculation of the EC effects on the fractal di- 
mensions T> q in the one and two-dimensional angular phase-space. We show that the self-similar 
structure of correlations is preserved beyond DLA and we give their analytic expression. As an 
application we calculate the EC effects on the fractal dimensions of QCD jets produced in e + —e~ 
collisions at LEP energy |7[ ^, [|. We find that these effects cannot be neglected either in the 
phenomenological analysis or in the theoretical considerations based on the perturbative QCD 
expansion. Under the EC effect, the fractal dimensions increase and, in the case of a running a s , 
strengthen the saturation phenomenon. In particular, the fractal dimension T>2 becomes now of 
order unity, in agreement with both numerical QCD simulations and experimental data. Note, 
from the theoretical point of view, the interesting interplay which appears with the property of 
KNO scaling of multiplicity distributions! 10 1. 

The plan of the paper is the following. In the nextcoming section 2, we recall the QCD 
formalism for the generating function of multiplicity moments for the global jet multiplicity, 
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formulate the EC problem and give an explicit solution for the first global moments. Section 
3 is devoted to the local multiplicity study of QCD jets, i.e. the correlations/fluctuations in 
small angular windows (angular intermittency[Q]). We compute the fractal dimensions and 
compare our analytical results (including rather significant EC effects) both with the relevant 
experimental data and with a computer simulation of QCD jets[|ll]]. Our conclusions can be 
found in section 4. 



2 The Global Jet multiplicity distribution 

We will use the following notations: Z(u, Q) stands for the QCD generating function of the 
global factorial multiplicity moments JF q (for a gluon of virtuality Q decaying into gluons). The 
total jet multiplicity (first global moment) is denoted N(Q). 

d q Z 

jF q =(„(„_!)..(„_ q + l)) jet = -g^\ u= v (6) 

When one takes into account energy conservation at the fragmentation vertex, Z is governed by 
an evolution equation which can be sketched as a classical fragmentation mechanism: 




where the black points represent the generating function of the jet and its sub-jets while the 
arrows stand for the jet and sub-jet axis directions. One obtains the well-known QCD evolution 
equation M : 



d\og(Q 



dZ(u,Q) f 1 2 dx . . 

~ / 7o — Z(u,Qx) Z(u,Q(l - x)) 
Jo x + 

= / ^—[Z(u,Qx)Z(u,Q(l-x))-Z(u,Q)] (7) 
Jo x 



where ^| + is the principal- value distribution coming from the triple gluon vertex and 6 is the 
angular aperture of the jet which plays the role of a time variable. The energy of the primordial 
parton is E and its virtuality is Q ~ Ed. 

Note that this equation can be enlarged to include other non-leading-log QCD contributions 
(the so-called Modified Leading-Log Approximation MLLA||). In fact an earlier study Q has 
shown that the EC effects play the crucial role in correcting the global moments. We shall thus 
focus our discussion on the EC effects but our method can be enlarged without major problems 
to the full MLLA equation. 

2.1 The mean multiplicity 

The solution for the mean multiplicity, N(Q) is known|| ||. let us rederive it in a language 
appropriate for further generalization to the correlation problem. Using Eq. (^), one obtains by 
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differentiation: 



dN(Q) 



2 dx 



Yo— (N(Qx) + N(Q(l - x)) - N(Q)) . 
x 



d\og{Q) 

Let us now distinguish the frozen-coupling regime and that with running a s . 

• Frozen coupling constant 

In this case, Eq. (||) can be solved using a power-like behaviour for N: 

N(Q) = iVoQ 7 

with the following link between 7 and the DLA value 70 of Eq.(||) : 



(8) 



7o 2 
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The function x is defined by: 

xO) 

, T , . <9iogr(z) 

where V(z) = ° — . When x 

oz 



1 " X(7) 
+ x) - 



(9) 



(10) 



0, X 



x 2 , 7 



70 — 7q/2 , which confirms that 



the EC correction to N is a (C(a s )) correction. However, this correction appears in 
the exponent of the multiplicity and is thus not negligeable at LEP energies. In practice, 
it " renormalizes" the multiplicity exponent, see Fig.|]. At LEP energy (E = Mz ), the 
correction can reach 20%. Note that x can De numerically replaced by x 2 up to x = 1.5. 
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Figure 1: a.70 as a function of the " renormalized" value 7; The grey band corresponds to the average 
value of 70 at LEP. 



Running coupling constant 

Let us now reconsider equation (^|) when a s is running. One has: 

a s (Q) cl 6 



7o = 4iV c 



2tt log (Q/A) 



; c 



^33 - 2n/' 
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where n/ is the number of flavors (n/ = 5 at LEP). The solution of Eq.(||) is of the form 
N(Q) = N exp (2<Vlog(Q)) , where c is a constant. Using the following general identity, 

N(Qf(x))/N(Q) ~ exp (2<yiog Q/(x)/A - 2 CN /logQ/A) 

~exp( C log/(x)/VIoiQ7A) = [f{x)] c l^ Q ' x . (11) 
one obtains from Eq.(||): 



N(Q) ~ iV exp (2<yiogQ/A) 

c 2 ^ C 2(l- x ( 7 )), (12) 

where 7 = c/ yTog Q/A. As already pointed out, x(t) — 7 2 when 7 is small, that is the 
EC effect on the exponent of the mean multiplicity is a non-leading correction. However, 
as in the frozen coupling case, it " renormalizes" the behaviour of the multiplicity. 



2.2 The global second moment 

If one goes to the second derivative of equation (0), one obtains the evolution equation for the 
global factorial moment, 1F 2 : 



dW 2 {Q) _ f 1 2 dx 

aiog(Q) 



f dx 

/ lo 2 —(lF 2 (Qx) + 2N(Qx)N(Q(l-x))+lF 2 (Q(l-x))-lF 2 (Q)). (13) 
Jo x 



The way of dealing with this equation is to assume that JF 2 (Q)/N 2 (Q) is a slowly varying 
function of Q. This property, which is one of the KNO scaling relations lF q (Q) oc N Q (Q), is 
known to be correct for the moments predicted by QCDQ. Using this structure together with 
the power-like behaviour of N, (see Eq. (S)) gives : 

«K^KCTFi*))) (Ar(<?>)2 - (14) 

The result is displayed on Fig.^. The EC effects are clearly seen as a serious damping of the 
KNO ratio 1F 2 /N 2 . The running of the coupling constant does not modify the result which is 
very close to the experimental dataQ (\bigf 2 in one hemisphere). One notices that Eq. ([l"4]) 
leads to a weak violation of KNO scaling due to the dependence of 7 on Q. 



3 The Local Jet multiplicity distribution 

Let us come now to the main topics of our paper, namely the computation of the local correlations 
between partons. The characteristic feature of the jet multiplicity structure is that the evolution 
equation for the local density of partons is linear. It can be deduced from the branching structure 
of jet fragmentation in QCD. 
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Figure 2: JF2 as a function of 7. The DLA result (4/3) is obtained when 7^0. The vertical grey band 
is for 70 at LEP , while the horizontal one corresponds to the experimental value quoted in ref. H] for 
one jet (in the forward hemisphere). 
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Let us denote by H (u, Q) the generating function for the density of factorial moments: 



d q H 



FJA) _ {n(n-l)..(n-q+l)) l 



(15) 



=1 A ~~ A 

The QCD evolution equation of the multiplicity in an observation window of size 6 pointing 
in a definite direction #0 from the jet axis is schematically described in the figure, there are 2 
contributions, one coming from the parton with a fraction x of the available energy, the other 
coming from the branch of the fragmentation process with 1 — x. Once substracted the variation 
of the phase-space from 6q to 9 by defining moments of the density as in Eq. (|l5[) , the only 
evolution comes from the energy degradation along the branch. In the DLA approximation, 
only one branch contributes to the hadron density evolution, since one parton keeps essentially 
undisturbed by the branching. In the EC case where we do not neglect the recoil effect, the two 
branches contribute. Note that one has to take into account the corresponding loss of virtuality 
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Qx and Q(l — x) respectively. One obtains : 
dH(u,Q) f 1 2 dx 



at 



f dx 

/ 7o 2 — (H(u, Qx) + H(u, Q(l - x))) 
Jo x + 

/ 7o 2 — ([H(u,Qx)-l] + [H(u,Q(l-x))-H(u,Q)]) (16) 
Jo x 



where log#o/$ = t plays the role of the evolution parameter and the phase-space splitting 
depends on the dimensionality, namely A oc 9 d . One branch of the iteration contributes 
to H(u,Qx) and the other to H(u,Q(l — x)) in the integrand, while the negative contribu- 
tions are necessary to ensure the finiteness and unitarity conditions, namely H(1,Q) = 1 and 
E A 9H(u, Q)/du\ u=l = E A *i(A) = N(Q). 

The solution of Eq.(|l6|) is obtained by using the KNO scaling property of the multiplicity 
distribution inside a cone of fixed angular aperture 6. In terms of the generating function H, it 
writes: 

H{u,XQ) = H{uN{QX)/N{Q),Q) (17) 
Inserting the property ([l?]) into Eq.(|l6|), and performing the q th derivative, one gets: 

- Hq i ( i — ^- i + ( - . " -' i - i 



Of, 



1 2 dx ff N(Qx) \ q ( N{Q{l-x)) 
70 x \\ N(Q) J + v. N(Q) 



H q f\ 2 -(x^+ (1 
Jo X 



\<rt 



(18) 



where we have used Eq. (|IT|) . Notice that Eq. (16), as well as (0), can be understood in terms of 
classical fragmentation models^]. The connection between the present formalism and classical 
fragmentation is also a consequence of the KNO relation (|i~7|). 

As a check of this procedure, let us recover the known results of the DLA approximation^]. 
Neglecting the EC terms in Eq.([l6|), one obtains : 



(19) 



dH(u,Q) f 1 2 dx 

jr. = / 7o — (H{u, Qx) - 1 , 

ot Jo X 



from wich one obtains the fractal dimensions 



3.1 Frozen coupling 

Equation (|D|) reads : 



1 dH q 
H q dt 



To 



(1 ~ x(Ql)) _ . 7 1 ~ x(Ql) 



(20) 



where one makes use of relation (pj)[). The corresponding corrections to the normalized factorial 
moments (|l|) and their fractal dimensions can be worked out : 



l + q + 



x(t) - x(<n) 



(21) 
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Figure 3: The multifractal dimensions 2?2,£>3,I>4 as functions of the coupling 7; The straight line is the 
DLA result for V4 as an example(||). The saturation limit is the dotted line at d = 1 . 

Notice that, while the correction to the unnormalized factorial moments, see Eq. (pl|), is a 
factor of order (1 — g 2 7 2 )/(l — 7 2 ) in the exponent and can be important, the normalized ones 
are less sensitive to EC corrections. Since x{l) — 7 2 up to 7 ~ 1.5, the fractal dimension (p0|) 
reads 

7(l + 7 2 )(g + l) ^ 1 - , 99 , 
V q ~ , q-y < 1.5, (22) 

leading to an enhancement factor with respect to the multiplicity exponent at the same order 
(o?J 2 )- However, when qj > 1.5, the modification becomes more important and q— dependent. 
One observes an important increase of the fractal dimension, see Fig.||. 

3.2 Running coupling 



Let us consider again Eqs.(|18|), ( po|) with 70 defined with the running coupling and 7 as in 
formula (|l2|). Neglecting both the variation of the coupling in the EC correction factors (higher 
order terms) and the x-dependence of the coupling (this approximation have been shown to be 
quite right in the range (z < .6) in ref. Q), one obtains : 



V q (z) ~ 1 



1 X(7) ~ X(g7) 

1-Q 



1 + VT 



(23) 



where 7 is expressed as in formula (|TC). Notice that the saturation effect is stronger due to the 
running of the coupling (see, e.g., Fig.| for the one-dimensional case). 

As a phenomenological application of our results, we perform a comparison of our predictions 
with experimental results, which while yet preliminary, have been presented as a first detailed 
analysis of correlations in a jet. As stressed in the experimental paper[13], the comparison of 



analytical QCD predictions with data on correlations is not easy due to the hadronization effects. 
The best candidate for a comparison is the study of particle correlations in a ring 9o±9 around 
the jet axis. The variation on 9 allows one to obtain the dependence in z, while considering as 
an observable the ratio log[J- q (z) / J- q (0)j minimizes the hadronization corrections^]. 

The results are displayed on Fig. ||. The result of a Monte-Carlo simulation of the QCD 



parton cascade is also shown on the same plot. The simulation is useful to take into account 
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Figure 4: log[jr g (z)/jT g (0)], as a function of the scaling variable, z = log[0o/#]/log[-E#o/A]. a: 
q = 3, with and without EC corrections; b: q = 2 (lower), q = 3 and q = 4 (upper), including 
EC corrections; the value 7 = .44 is consistent with the LEP value of a s . 



the corrections with respect to analytical calculations which occur when the parton cascade 
is correctly reproduced, e.g. overlap effects between angular domains, subasymptotic energy- 
momentum effects, hadronization effects (in part) etc.... For this comparison we make use of 
the following parameters : a s (E6o) = .135 for #0 = 25° (corresponding to a s = .12 for the whole 
jet with Qq = 45 Gev). This leads to 7 = .44 or, from Eq.(10), 70 = .5. 

The analytical QCD predictions happen to be in reasonable agreement with both experi- 
mental results and numerical Monte-Carlo simulations, taken into account the approximations 
considered in the calculation. In particular, some bending of the factorial moments at small 
aperture angle 9 is qualitatively reproduced. This effect is a typical predictions of the running 
coupling scheme, which is thus favoured by the data with respect to the frozen coupling case. 

The full study of hadronization effects remain to be done. However, we may notice that 
when the slope of the intermittency indices (Fig. ||]) goes to zero, the effective coupling 7 is 
around 1 and one enters the not well known strong coupling regime of QCD. As a matter of 
fact, the only free parameter left in the Monte-Carlo simulation is the effective cut-off of the 
parton cascade i.e. the maximum value a s , oto, one considers in the development of the parton 
branching process. In the example considered in the figure, we have found ao — -54, which is a 
quite reasonable value for switching on hadronization. All in all it is a positive fact that Local 
Parton-Hadron Duality seem to work reasonably for such refined quantities as correlations, once 
some caution is taken to minimize non-perturbative QCD corrections. 



4 Conclusions 



The Energy Momentum Conservation effects (or recoil effects) on the correlation properties 
of QCD jets are thus under some analytical control. The method proposed to handle these 
effects is the search for scaling solutions of the QCD generating function of local multiplicity 
density moments in angular phase-space. It is interesting to note the deep connection between 
the KNO-scaling[10] and the self-similar properties of correlations in the jet. This structure is 
reminiscent of general properties of fragmentation models and could be more fundamental than 
the level of approximation used to derive the result from QCD. Indeed, the same qualitative 
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Figure 5: Phenomenological predictions for log[F q (z) /F g (0)] as a function of the scaling variable z. 
Upper curves, triangles and circles: q — 3, lower ones q = 2. Continuous lines: the analytical calculations 
from [|; circles: Preliminary Delphi data from Ref . |]i~3f ; squares: Numerical data from the Monte-Carlo 
simulation fi~l|| 



result is expected from the solution of the equation including EC terms plus other dominant 
next-to-leading effects (e.g. MLLA approximation scheme, see ref.(||)). 

The phenomenological outcome of this study is that the analytical predictions come close 
to the experimental data and confirmed by the numerical simulation. Parton-Hadron dual- 
ity for correlations seems to be supported provided hadronization effects can be minimized. 
Thus the self-similar structure of the QCD branching processes may eventually emerge from the 
background in the experimental analysis. Hadronization becomes dominant when the factorial 
moments bend and acquire a zero slope (z ~ 0.5 in terms of the scaling variable of formula (g)). 
If one would insist to increase z by reducing the opening angle of the detection window, one 
would find properties closer to soft processes, and thus difficult to predict theoretically. 

From the theoretical point of view, the partial cancellation of EC corrections observed in the 
scaled factorial moments, leading to a moderate correction to the behaviour determined by the 
leading-log approximation, is a consistency check of the perturbative resummation predictions. 
Note however that this remark does not hold for moments of high rank q > 4, where a stronger 
correction is expected. The question remains to know whether this indicates a limitation of the 
perturbative calculations or the signal of a dynamical mechanism. 
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